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Abstract 


Three light sterile neutrinos (vé, v} and v3) are introduced to accommo- 
date all the available neutrino data: the atmospheric neutrino anomaly is 
explained by v, — vj, oscillation with maximal mixing; the solar one is due 
to ve — vé oscillation of small angle Mikheyev-Smirnov- Wolfenstein type; the 
Liquid Scintillation Neutrino Detector data is from ve — vy, oscillations, so 
that the neutrinos can be the hot component of the dark matter. The big 
bang nucleosynthesis constraint is satisfied by taking the tau neutrino to be 
10 MeV heavy. The v, decay is discussed in a model of gauge mediated super- 
symmetry breaking. The decay mode v, > Gy with G being the gravitino is 
proposed. The v, has a rather long lifetime ~ 10° — 10!° sec. Its implication 


to the Gamma-ray Burst is discussed. 
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I. INTRODUCTION 


Various experiments have provided growing evidence that neutrinos are massive. First 
the recent data on the atmospheric neutrinos from Super-Kamiokande experiment has shown 
neutrino oscillations |I]. It implies that the p-type neutrino has a maximal mixing with 
other neutrino x (x # e), and the mass squared difference is Am2, ~ 3 x 107? eV’. The 
solar neutrino deficit problem can be explained by either the Mikheyev-Smirnov- Wolfenstein 
(MSW) solution or the vacuum oscillation BJ. The MSW solution allows two sets of 
parameters: Am, ~ 5 x 10~® eV? with sin?26., ~ 6 x 107-3, and Am2, ~ 2 x 107° 
eV? with sin? 20e, œ 0.8. The vacuum oscillation solution is Am2, ~ 8 x 107" eV? with 
sin? 20a) ~ 0.8 H If the direct observation of D, — De oscillations at the Liquid Scintillation 
Neutrino Detector (LSND) experiment |H] is considered, relevant parameters are Am, ~1 
eV? and sin?(20.,,) ~ 1072. 

In addition, astrophysics and cosmology also give us some information for neutrino 
masses. Various measurements support the existence of the dark matter. One scenario 
for the dark matter is that not all of the cosmological dark matter is cold, there is some 
hot component [6]. Canonical candidate for the hot dark matter (HDM) is the neutrino 
with mass at the scale of electron Volt |]. If massive neutrinos are stable, their masses are 
limited to be less than few tens of eV from the standard cosmology, in order to avoid the 
over-closure of the Universe [f]. 

Non-vanishing neutrino mass might be the first discovery of the physics beyond the stan- 
dard model (SM). The mass pattern of neutrinos is valuable information for the exploration 
of the physics related to the flavor puzzle in the SM. Theoretically, several mechanisms have 
been suggested to accommodate massive neutrinos. One of the most popular scenarios is the 
seesaw mechanism which naturally explains the smallness of neutrino mass by introducing 
heavy right-handed neutrinos [B]. Another example is the supersymmetric extension of the 
SM with R-parity violation. The trilinear lepton number violating interactions induce small 
neutrino masses at loop level JJ. However, the large mixing of neutrinos is not expected 
from our experiences of the other SM fermions. 

Phenomenological analyses suggest the existence of light sterile neutrinos. As least one 
light sterile neutrino is necessary to explain all the three different scales of Am? [I0]. It is 
interesting and more natural that there are three light sterile neutrinos v? (i = e, y, T). Then 
the first two generations can fully explain all the above neutrino data. In this framework, the 
neutrino mixing pattern is uniquely fixed. The atmospheric and solar neutrino anomalies 
are due to the v, — v} and ve — vå oscillations, respectively. The large mixing between v, 
and vi can be understood naturally. And the LSND evidence is attributed to the ve — v, 
mixing. At least one of the neutrino pairs should have mass around few eV, implying that 
neutrinos can be the hot component of the dark matter. This scenario was studied before 
E. In this paper we reconsider it in a different theoretical background. 

In this scheme, more physics can be addressed. The third generation neutrinos have a 
separate story. The above requirement for the first two generations leaves a wider room 


‘Recently the low mass solution has been revived which is Am?, ~ 1077 eV? with large mixing 


g. a 


for the v,. The v, can be either ordinary with a mass smaller than 10 eV, or exotic, for 
instance, with m,, ~ (1 — 10) MeV [Ll]. The physics of tau-neutrinos is influenced by 
the astrophysical constraint from the big bang nucleosynthesis (BBN): the number of light 
neutrino species in thermal equilibrium at the BBN era is limited to N, < 4.2 [I3]. It should 
be seriously examined in a scenario with more than four light neutrinos. It is to be noted 
that the small mixing angle for the ve — vå oscillations causes only the active neutrinos to 
contribute at the BBN era [f[4]. The oscillation time is too long to make the vå in thermal 
equilibrium. In fact, the analyses for the solar neutrino data have shown that the only viable 
pattern is the small angle MSW type [I5]. Therefore the first two generations contribute 
a factor three to N,. As for the v, and v$, at least one of them should decouple at the 
time of BBN. If v, is light (my, < 10 eV), the sterile tau neutrinos must decouple from 
other particles in the primeval plasma. This is equivalent to the case of introducing just 
two sterile neutrinos [f[G]. On the other hand, we may consider m,, > 1 MeV so that v, 
decouples [fT]. We are interested in this latter situation which has alternative theoretical 
motivation {7}. In addition, the economy of the first two generation oscillation scenario 
becomes more meaningful. 

It is a theoretically challenging problem to construct a spectrum of very light sterile 
neutrinos. Several ideas have been suggested. They are related to extra symmetry [8], 
mirror world [9], axino [20], modulino [Pi], gravitino [PJ], string theory PJJ, gauge medi- 
ated supersymmetry breaking (GMSB) [24], U(2) symmetry [25], seesaw mechanism P6), 
compositeness 27], top-flavor model and extra dimensions [29]. 

This paper is organized as follows. In the next section, we present a mass matrix which 
is consistent with all the available neutrino data. The naturalness of the assumption is 
discussed. In Sect. III, the decay of the tau neutrino is studied in a supersymmetric model. 
The decay mode v, — Gy is proposed. Further discussions including the gamma-ray burst 
(GRB) are made in Sect. IV. The summary is given in the last section. 


Il. NEUTRINO MASS MATRIX 


We assume that there are three light sterile neutrinos. For each generation of fermions, 
there is an associated sterile neutrino. As will be seen in the following, the atmospheric 
and LSND data will be accommodated naturally. And the tau neutrino can be as heavy 
as (1 — 10) MeV without contradicting experiments. The full 6 x 6 neutrino mass matrix 
M consists of three parts. They are 3 x 3 active Majorana neutrino mass matrix M2ve, 
3 x 3 Dirac mass matrix MP1 which stands for the mixing between the active and sterile 
neutrinos, and 3 x 3 Majorana mass matrix of the sterile neutrinos, respectively, 


Mactive Mbirac 
M = (ype on) : (1) 


For the M?cv*, we assume that except for the v,v, component it originates from the 
ordinary seesaw mechanism [$]; and the v,v, component is about (1—10) MeV. The ordinary 
seesaw mechanism suppresses neutrino masses to be v?/Mgyr, where v is the electroweak 
scale and Maur the grand unification (GUT) scale so that v?/Mcaur ~ 107° eV. This is 
natural because there are certain evidences that the new physics beyond the SM is the GUT 


BO]. The introduction of the sterile neutrinos does not affect the gauge coupling unification. 
Furthermore, as can be seen below, this 107? eV mass just fits the atmospheric neutrino 
data. In other words, the neutrino mass itself is also indirect evidence of the GUT. Note 
that even without the heavy right-handed neutrinos, the seesaw mechanism still works in 
the GUT scenario. It is from the general argument based on the effective field theoretical 
point of view. The MeV tau-neutrino has been discussed in Ref. |I]. We introduce it 
here to avoid the BBN constraint and to make the two generation oscillation scenario more 
economical. 

For the M?*°, we adopt that it is proportional to the up-type quark masses. This simi- 
larity is understandable from the GUT point of view, if the sterile neutrinos are regarded as 
right-handed neutrinos. It is remarkable that this simple proportionality can be compatible 
with all the current neutrino data. In this quark mass matrix, the up quark is taken to be 
10-4m,. The up quark mass has not yet been determined experimentally. Even a massless 
up quark is allowed [B] BJ]. 

For the M%teile we simply assume that it is approximately a vanishing matrix. This 
assumption makes the point of light sterile neutrinos. It could be relaxed to certain extent 
that the Merile plays no other roles in understanding the neutrino experiments. 

Since the main understanding of the neutrinos does not involve intergenerational flavor 
mixing, without losing generality, we further assume that all these mass matrices are diag- 
onal. We make a remark here that the other zero entries of the Majorana mass matrices 
should be regarded as ~ 1075 eV in general. This is inevitable from the consideration of 
quantum gravity B3]. However, such a number is not essential for the understanding of the 
neutrinos in this framework. 

Therefore the neutrino mass matrix M takes the following form, 


y2 
O (=) 0 0 AM, 0 0 
v2 
0 O ) 0 Am. O 
M= Meur (2) 
0 0 10 MeV 0 0 AM 
AMy 0 0 0 0 0 
0 AMe 0 0 0 0 
0 0 AM: 0 0 0 


where A is a constant. 

Here starts numerical analyses. As has been mentioned, v? / Mayur is taken to be 3 x 1073 
eV. We take Am, œ 1 eV. It is easy to see that v, and Vi have an almost maximal mixing with 
Am? ~ 3 x 107 eV’. This explains the atmospheric neutrino data. The mixing between 
Ve and vê is very small: sin? 20 ~ [2 x 1074/(3 x 1073)|? ~ 4 x 1073 with Am? ~ 107° eV?. 
This is the small angle MSW solution for the solar neutrino data. The Dirac mass of tau 
neutrino is 100 eV. This leads to no observable physical consequences of sterile tau neutrinos 
because the Majorana mass of v, is about 10 MeV which makes the v, — v mixing to be 
at ~ O(10~°). As discussed before, this small mixing precludes the sterile tau neutrinos 
from contributing at the BBN era. The LSND data of the ve — v, oscillations is accounted 
by Am?2,, ~ 1 eV’. The intergenerational mixings of the neutrinos are determined by the 
neutrino mass matrix and the charged lepton mass matrix, which is too complicated to be 
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fixed in this framework. Since the Dirac mass of the electron neutrino is negligible, the 
Ve — V, mixing depends only on M**’* and the charged lepton mass matrix. In the case 
where M**ve is diagonal in certain weak basis, we expect that the charged lepton mass 
matrix will fix the Ve — v, mixing, namely m./m, ~ 5 x 107? which is consistent with the 
data. The v, and vf have masses around 1 eV so that they can be the HDM. The v, cannot 


H 
be stable due to cosmological considerations. 


II. THE DECAY OF pv, 


In the previous section, phenomenology of three light sterile neutrinos is discussed, and 
a heavy Majorana tau neutrino is necessarily assumed. To study the 10 MeV tau neutrino 
and its decay, a theoretical model needs to be specified. First, this v, may decay via the 
ordinary weak interaction, namely into ete ve through the W-boson exchange. Its lifetime 
is 
19273 1 


xv —— v 02.3 x 1—5. 3 
e oe e ~ °° * Gm, MN e * 8) 


Ty 
Since the neutrino mass matrix leads to a very small mixing, only the charged lepton mass 
matrix will determine V,,. It is generally expected that |V.,|? ~ Mme/Mm, ~ 3 x 1074 which 
results in 7,, ~ 10° sec. However this decay mode is disfavored by the observation of the 
Supernova 1987A 4]. Even if this decay is marginally viable by assuming smaller Ver, 
something more about the v, decay can be addressed. 

There is an alternative decay channel which is irrelevant to the e — r mixing. It is 
inherent in the model of the (1 — 10) MeV v,. In Ref. [B5], a supersymmetric model for the 
charged lepton masses is suggested. In this model, because of family symmetries, 7 lepton 
gets mass from the nonvanishing vacuum expectation value (VEV) of the Higgs field. The 
muon gets mass only due to the sneutrino VEV which violates the family symmetries. The 
lepton number violation is allowed in the model. The sneutrino VEV, which generates a 10 
MeV mass to the tau-neutrino, is as large as (5 — 10) GeV, resulting in a mixing between 
the ordinary neutralinos and the v, [l7]. This model adopts the framework of the GMSB 
Boj. The v, in fact is the lightest neutralino (except for the gravitino) and the 7 the lightest 
chargino. Therefore v, can decay to Gy where the gravitino has mass of few eV’s. From 
the neutralino mass matrix in Ref. [I7], it can be seen that v, mixes with Bino at the level 
of m,y/M; ~ 107? where m,y is the mixing mass of the neutrino and neutralinos. In 
ordinary GMSB models BE] the next-to-lightest neutralino, photino 7, decays to G and y. 
The relevant matrix element (G'y|7) can be written in the same way as in current algebra 


X gs 1 zyja 
(Gyl) = plal) 5 (4) 
where the supersymmetry current is 
jh = PYG + oF By + (5) 


with VF being the supersymmetry breaking scale ~ 100 TeV. From the second term in 
the above equation, the matrix element can be easily evaluated. In our model, through 
considering the suppression factor (m,4/M;)°, the v, + Gy decay rate is then given by 
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This is equivalent to a lifetime of 10!° sec, if we switch off the V.,. Even though this lifetime 
is very long compared to the BBN epoch, it is still short enough compared to the age of the 
Universe. 

The BBN constraint for the number of light neutrino species is counted by ve, V, and v$. 
A lifetime of 10° sec for v, is at the finishing stage of the BBN. Such tau-neutrinos may be 
considered as unstable in the BBN time. With ve being the decay product, it is probably 
helpful for the understanding of the Deuterium production and the structure formation in 
the cosmology [B7]. It is however more appropriate to take v, with lifetime of > 10° sec as 
stable in the BBN epoch B9]. In this case, v, decouples in the BBN time. And it must be 
Majorana neutrino B7B8} as we have assumed. 


IV. DISCUSSIONS 


Several discussions are necessary about our scenario for neutrinos. The essential new 
ingredients of this scenario are the 10 MeV Majorana tau-neutrinos and the three light 
sterile neutrinos. A theoretical model of the 10 MeV v, and its decay have been discussed. 
However, the lightness of the sterile neutrinos is still lack of theoretical discussion. To our 
knowledge, for more than one light sterile neutrinos, there exists two attractive frameworks. 
One is the compositeness idea [27]. And the other comes from the mirror world [19]. They 
are worthy to be studied further. We introduce them phenomenologically in this paper. 

While the R-parity violating supersymmetry is used to generate the MeV mass for v, 
it would be expected that the trilinear R-parity violating interactions contribute to the 
other entries of M?*“v°, However, this contribution is rather small. It arises at one-loop 
level, where a soft term of the relevant trilinear R-parity violating interaction, A-term, is 
involved. In the GMSB, A-terms are very small. The contribution is in general smaller than 
1073 eV. 

Experimentally, our scenario can be tested in the very near future. First, the current 
laboratory mass limit for v, is 18.2 MeV |B9]. It can be improved to be a few MeV in the 
near future, which could be a serious challenge to this scenario. 

Second, future neutrino experiments will determine whether the atmospheric neutrino 
deficit is due to v, — vj, or V, — Vr oscillation. Although the present data favor the v, — v, 
oscillation, the v, — v} oscillation cannot be ruled out on the basis of the global fit to the 
full set of observable 0]. Recent Super-Kamiokande experiment claims that the v, > vi 
oscillation for the atmospheric neutrino anomaly is ruled out at 99% C.L. Hl]. However, it 
has been argued that more data and analyses are necessary before the conclusion is made 
EJ. The long baseline neutrino experiments [#3] will check the oscillation channel. On the 
other hand, the oscillation channel for the solar neutrino deficit can be hopefully picked out 
by the new SNO observation M4] together with previous experiments. 

Third, the short baseline neutrino experiment BooNE will confirm whether LSND 
anomaly is real. Even though the KARMEN [45] group has reported that a large part 
of the favored parameter region of LSND is excluded, full confirmation of the LSND results 
still awaits future experiments. For example, some of the LSND solution space is preserved 


by a combined statistical analysis of both the LSND evidence and the KARMEN exclusion 
HG. 

Cosmological and astrophysical aspects of this scenario deserve further studies. The 
implications on the BBN and HDM have been discussed in the previous sections. For the 
decay mode v, — G4, it might be the partial reason of the GRB EJ. In the process of the 
stellar collapse, a huge number of tau neutrinos are emitted. The decay product y is of few 
MeV or 100 keV which is the typical energy scale of y in the GRB. Suppose that the lifetime 
of v, is of order 1018 sec ( the Ve, is negligible small < 1077), 10 MeV v, emitted from the 
Supernova 1987A did not have enough time to decay when they arrived in the Earth #8). 
When the supernova occurs at cosmological distance (> 1076 — 10?” cm) from us, however, 
the GRB will be detected. The (1 — 10) MeV v, may help our understanding of the cold 
dark matter (CDM) [9]. As for the sterile neutrinos, it will be harmless for our numerical 
results to assume that vê has Majorana mass of thousand electron Volts. Then this v° is 
also the candidate of the warm dark matter (WDM) [0]. 

On the other hand, if v, purely decays to Gy with a lifetime ~ 10!° sec, too much 
cosmic background radiation and matter density are caused ff]. For reducing the background 
radiation, one way maybe to prolong the v, lifetime from 10!° sec until to 10? sec. This can 
be achieved through adjusting the supersymmetry breaking scale (along with some other 
model parameters) up to 104 TeV. (In this case, the gravitino mass rises from few eV to 
100 keV.) To avoid the matter density problem, we may need some non-standard cosmology. 
For example, in a recent study [61], it has been shown that if the largest temperature of the 
radiation era is as low as 0.7 MeV, even a stable 10 MeV tau neutrino is consistent with the 
cosmology. 

Within the framework of GMSB, the superpartners of the sterile neutrinos can be very 
light because they are not involved in the gauge interactions. However, their mixings with 
the ordinary sneutrinos which are typically 100 GeV — 1 TeV heavy, are very small. So their 
lightness has no influence on the physics of the BBN. Meanwhile, these light scalar particles 
can be also the candidate of the CDM pJ. 

Compared to the previous studies of MeV tau neutrinos, the distinguished feature of 
this scenario is that the (1 — 10) MeV tau neutrino is Majorana particle and its decay 
products do not contain Majoron. Different from the pseudo-Dirac neutrino scenarios [12], 
our explanation of the solar neutrino deficit is due to small angle MSW solution. The 
previous motivation was mainly to reduce the number of neutrino species in BBN era to be 
smaller than 3, the present one is to keep that number just as 3. 


V. SUMMARY 


Before summarizing the main points, let us list the dark matter candidates in this model. 
Besides the v, and v$ which have been taken as the HDM, there are many candidates for the 
CDM and the WDM. They are (1—10) MeV v,’s, keV massive v/°’s, light sterile sneutrinos 1? 
and gravitinos. However, these are just possibilities. Because the lifetime of v, has not been 
really fixed which still ranges from (10? — 1013) sec (or longer); the heaviness of vê and 7 are 
just assumed, which need to be studied from the theoretical mechanism for the lightness of 
sterile neutrinos. Note that by taking the supersymmetry breaking scale as ~ 100 TeV, the 
gravitino is not the dark matter in simple GMSB models. However by raising the VF, keV 


G can be the WDM (regardless the R-parity violation). In addition, in the GMSB scenario, 
lightest stable baryons with mass of ~ 100 TeV in the hidden sector is also a candidate 
for the CDM [6]. Among the five possibilities, the long-lived and decaying v, which is 
specific in this model, is reasonably favored to be the CDM in this model according to our 
discussion. For the stable v,, the matter density was discussed in Ref. [Bil]. It is expected 
that the same discussion can be made for the very long-lived v,. The requirement Q,,h? < 1 
implies m,. > 3 MeV for the reheating temperature being 0.7 MeV. 

In summary, we have introduced three light sterile neutrinos v3, vj, and v$ to accom- 
modate all the available neutrino data. The atmospheric neutrino anomaly is explained by 
Vu — Va oscillation with maximal mixing and the solar one by ve — v$ oscillation of small 
angle MSW type. The LSND data of the ve — vy, oscillations with Ae, ~ 1 eV? can be ac- 
commodated, as providing canonical candidate of the HDM. The BBN constraint is satisfied 
by taking the tau neutrino to be 10 MeV heavy. The v, decay has been discussed in a model 
of GMSB. The decay mode v, + Gy has been proposed. The implication to the GRB has 
been pointed out. Although the current data do not favor this oscillation scenario, it will be 
conclusively checked by atmospheric neutrino experiments, long baseline experiments and 
solar experiments in the near future. 


S 


ACKNOWLEDGMENTS 


We would like to thank Carlo Giunti, Tan Lu for helpful discussions, and Xin-min Zhang 
for helpful discussions, drawing our attention to Ref. and informing us that a similar 
scenario in the mirror world model with left-right symmetry was considered by R. N. Mo- 
hapatra and X.-M. Zhang. This work was supported in part by the Fund of High Energy 
Physics Development of China, BK21 Program of the Ministry of Education of Korea and 
the National Natural Science Foundation of China. 


REFERENCES 


Y. Fukuda et al. (Super-Kamiokande Collaboration), Phys. Rev. Lett. 81 (1998) 1562. 
S. P. Mikheyev and A. Yu Smirnov, Sov. J. Nucl. Phys. 42 (1985) 913; L. Wolfenstein, 
Phys. Rev. D 17 (1978) 2369. 

3| For an example, see V. Barger, R. J. N. Phillips and K. Whisnant, Phys. Rev. Lett. 69 

(1992) 3135. 

4| For a recent review, see M. C. Gonzales-Garcia and C. Pena-Garay, |hep-ph/000904] 

5] C. Athanassopoulos et al. (LSND Collaboration), Phys. Rev. Lett. 75 (1995) 2650. 

6] J. R. Primack, J. Holtzman, A. Klypin, and D. O. Caldwell, Phy. Rev. Lett. 74 (1995) 

2160; For a review, see J. R. Primack and M. A. K. Gross, fAstro-ph/0007165]. 

7| For a review, see E. W. Kolb and M. S. Turner, The Early Universe (Addison-Wesley, 

Reading, 1990). 

8] M. Gell-Mann, P. Ramond and R. Slansky, in Sanibel talk, CALT-68-709 (Feb. 1979), 
and in Supergravity (North Holland, Amsterdam, 1979); T. Yanagida, in Proc. of the 
Workshop on Unified Theory and Baryon Number of the Universe (KEK, Tsukuba, 
1979); R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44 (1980) 912. 

[9] C. Aulakh and R. Mohapatra, Phys. Lett. B 119 (1982) 136; S. Dimpoulos and L. Hall, 
Phys. Lett B 207 (1988) 210; C. Liu, Mod. Phys. Lett. A 12 (1997) 329; M. Drees, S. 
Pakvasa, X. Tata and T. Veldhuis, Phys. Rev. D 57 (1998) R5335; E. J. Chun, S. K. 
Kang, C. W. Kim and U. W. Lee, Nucl. Phys. B544 (1999) 89; G. Bhattacharyya, H. V. 
Klapdor-Kleingrothaus and H. Pas, Phys. Lett. B 463 (1999) 77; M. Bisset et al., Phys. 
Rev. D 62 (2000) 035001; A. Abada and G. Bhattacharyya, Phys. Rev. D 63 (2001) 
017701. 

[10] D. O. Caldwell and R. N. Mohapatra, Phys. Rev. D 48 (1993) 3269; J. T. Peltoniemi 
and J. W. F. Valle, Nucl. Phys. B 406 (1993) 409; J. J. Gomez-Cadenas and M. C. 
Gonzales-Garcia, Z. Phys. C 71 (1996) 443; N. Okada and O. Yasuda, Int. J. Mod. 
Phys. A 12 (1997) 3669; S. M. Bilenky, C. Giunti and W. Grimus, Eur. Phys. J. C 1 
(1998) 247; Q. Y. Liu and A. Yu Smirnov, Nucl Phys. B 524 (1998) 505; V. Barger, S. 
Pakvasa, T. J. Weiler and K. Whisnant, Phys. Rev. D 58 (1998) 093016; E. J. Chun, C. 
W. Kim and U. W. Lee, Phys. Rev. D 58 (1998) 093003; Y. Chikira, N. Haba and Y. 
Mimura, fhep-ph/9808254; Q. Shafi and Z. Tavartkiladze, Phys. Lett. B 448 (1999) 46; 
V. Barger, Y.-B. Dai, K. Whisnant and B.-L. Young, Phys. Rev. D 59 (1999) 113010; 
C. Liu and J. Song, Phys. Rev. D 60 (1999) 036002. 

[11] N. Okada, Phys. Lett. B 388 (1996) 106; D. Ring, Phys. Rev. D 55 (1997) 5767; 

X.-M. Zhang, talk at the Annual Meeting of High Energy Physics of China (Chengdu, 
1999). 

[12] M. Kobayashi, C. S. Lim and M. M. Nojiri, Phys. Rev. Lett. 67 (1991) 1685; C. Giunti, 
C. W. Kim, and U. W. Lee, Phys. Rev. D 46 (1992) 3034; A. Geiser, Phys. Lett. B 444 
(1998) 358; D. Chang and O. C. W. Kong, Phys. Lett. B 477 (2000) 416. 

[13] For a review, see K. A. Olive, Lectures at the Advanced School on 
Cosmology and Particle Physics (Peniscola, June 1998). 

[14] R. Barbieri and A. Doldov, Nucl. Phys. B 429 (1991) 743. 

[15] For recent analyses, see J. N. Bahcall, P. I. Krastev and A. Yu Smirnov, Phys. Rev. D 

58 (1998) 096016; M. C. Gonzales-Garcia, P. C. de Holanda, C. Pefia-Garay and J. W. 

F. Valle, Nucl. Phys. B 573 (2000) 3. 


V = 


16 


17 
18 
19 


20 


21 
22 
23 
24 
25 
26 
27 
28 
29 


30 


31 


32 


33 


34 


35 
36 


37 


W. Królikowski, Acta Phys. Polon. B 30 (1999) 227; D. Suematsu, Int. J. Mod. Phys. 
A 15 (2000) 3967. 

C. Liu and H. S. Song, Nucl. Phys. B 545 (1999) 183. 

E. Ma and P. Roy, Phys. Rev. D 52 (1995) R4780. 

R. Foot and R. R. Volkas, Phys. Rev. D 52 (1995) 6595; Z. G. Berezhiani and R. N. 
Mohapatra, Phys. Rev. D 52 (1995) 6607. 

E. J. Chun, A. S. Joshipura and A. Yu. Smirnov, Phys. Rev. D 54 (1996) 4654; E. J. 
Chun, Phys. Lett. B 454 (1999) 304; 

K. Benakli and A. Yu Smirnov, Phys. Rev. Lett. 79 (1997) 4314. 

T. Lee, hep-ph/9902396 

P. Langacker and J. Wang, Phys. Rev. D 58 (1998) 115010. 

G. Dvali and Y. Nir, JHEP 10 (1998) 014. 

L. J. Hall and N. Weiner, Phys. Rev. D 60 (1999) 033005. 

B. Brahmachari, Phys. Lett. B 461 (1999) 243. 

N. Arkani-Hamed and Y. Grossman, Phys. Lett. B 459 (1999) 179. 

E. Malkawi, E. I. Lashin and H. Widyan, Phys. Rev. D 62 (2000) 033005. 

R. N. Mohapatra and A. Perez-Lorenzana, Nucl. Phys. B 576 (2000) 466; Y. Grossman 
and M. Neubert, Phys. Lett. B 474 (2000) 361 and references therein. 

P. Langacker and M.-X. Luo, Phys. Rev. D 44 (1991) 817; C. Giunti, C. W. Kim and U. 
W. Lee, Mod. Phys. Lett. A 6 (1991) 1745; U. Amaldi, W. de Boer and H. Furstenau, 
Phys. Lett. B 260 (1991) 447; J. Ellis, S. Kelly and D. Nanopolous, Phys. Lett. B 260 
(1991) 131. 

H. Georgi and I. N. McArthur, HUTP-81/A011 (1981); 

K. Choi, C. W. Kim and W. K. Sze, Phys. Rev. Lett. 61 (1988) 794. 

C. Liu, Int. J. Mod. Phys. A 11 (1996) 4307; Liu and Song, in [0]; J. E. Kim, B. Kyae 
and J. S. Lee, Phys. Lett. B 477 (1999) 110. 

E. Kh. Akhmedov, Z. G. Berezhiani, G. Senjanovic and Z.-J. Tao, Phys. Rev. D 47 
(1993) 3245; For reviews, see E. Kh. Akhmedov, Lectures at Trieste 
Summer School in Particle Physics (June, 1999); E. Witten, hep-ph/0006333, 

A. Dar and S. Dado, Phys. Rev. Lett. 59 (1987) 2368; R. N. Mohapatra, S. Nussinov 
and X. Zhang, Phys. Rev. D 49 (1994) 3434. 

D. Du and C. Liu, Mod. Phys. Lett. A 8 (1993) 2271; A 10 (1995) 1837. 

For reviews, see G. F. Giudice and R. Rattazzi, Phys. Rep. 322 (1999) 419; S. L. 
Dubovsky, D. S. Gorbunov and S.V. Troitsky, fhep-ph/9905466. 

For a review, see G. Gyuk and M. S. Turner, Nucl. Phys. B (Proc. Suppl.) 38 (1995) 
13; For recent studies, see A. D. Dolgov, S. H. Hansen, S. Pastor and D. V. Semikoz, 
Nucl. Phys. B 548 (1999) 385; M. Kawasaki, K. Kohri and K. Sato, Phys. Lett. B 430 
(1998) 132. 


[38] For a review, see G. G. Raffelt, Nucl. Phys. B (Proc. Suppl.) 70 (1999) 169; For a recent 


study, see J. B. Rehm, G. G. Raffelt and A. Weiss, Astron. Astrophys. 327 (1997) 443 
and references therein. 


[39] D. E. Groom et al. (Particle Data Group), Eur. Phys. J. C 15 (2000) 1. 
[40] For a recent analysis, see N. Fornengo, M. C. Gonzales-Garcia and J. W. F. Valle, Nucl. 


Phys. B 580 (2000) 58. 


[41] S. Fukuda et al. (Super-Kamiokande Collaboration), Phys. Rev. Lett. 85 (2000) 3999. 


10 


[42] R. Foot, fhep-ph/000700% 
M. C. Gonzales-Garcia, talk at ICHEP 2000 (Osaka, August, 2000). 


[43] C. Yanagisawa (KEK-SK Collaboration), Proc. Int. Workshop on Phys. beyond the 
Standard Model: from Theory to Experiment (Valencia, 1997), eds. I. Antoniadis, L. 
Ibanez and J. W. F. Valle (World Scientific, 1998); 

D. Michael et al. (MINOS Collaboration), Proc. XVI Int. Workshop on Weak Interac- 
tions and Neutrinos (Capri, 1997), eds. G. Fiorilo et al. [Nucl. Phys. B (Proc. Suppl.) 
66 (1998) 432]; 

T.-H. Ho, talk at the ITP Academic Meeting (Beijing, 2000). 

44] SNO Collaboration, Phys. Lett. B 194 (1987) 321; S. M. Bilenky and C. Giunti, Phys. 

Lett. B 320 (1994) 323. 

45] B. Zeitnits, talk at Neutrino ’98 (Takayama, June 1998). 

46| K. Eitel, New Jour. Phys. 2 (2000) 1. 

47] For reviews, see S. Blinnikov, T. Lu, fstro-ph/0002178. 

48] A. H. Jaffe and M. S. Turner, Phys. Rev. D 55 (1997) 7951. 

49] S. Dodelson, G. Gyuk and M. S. Turner, Phys. Rev. Lett. 72 (1994) 3754. 


50] M. Patel and G. M. Fuller, hep-ph/0003034. 


51] G. F. Giudice, E. W. Kolb and A. Riotto, hep-ph/0005123} 
52] I. I. Tkachev, Sov. Astron. Lett. 12 (1986) 315. 


11 


